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The scale-free (SF) structure that commonly appears in many complex networks is one of the 
hot topics related to social, biological, and information sciences. The self-organized generation 
mechanisms are expected to be useful for efficient communication or robust connectivity in socio- 
technological infrastructures. This paper is the first review of geographical SF network models. We 
discuss the essential generation mechanisms to induce the structure with power-law behavior and 
the properties of planarity and link length. The distributed designs of geographical SF networks 
without the crossing and long-range links that cause the interference and dissipation problems are 
very important for many applications such as in the Internet, power-grid, mobile, and sensor systems. 

PACS numbers: 89.75.Fb, 89. 75. Da 



I. INTRODUCTION 

As a breakthrough in network science [1], it has been 
found [2] that many real systems of social, technolog- 
ical, and biological origins have the surprisingly com- 
mon topological structure called small-world (SW) [3] 
and scale-free (SF) [4]. The structure is characterized 
by the SW properties that the average path length over 
all nodes (vertices) is short as similar to that in ran- 
dom graphs, and that the clustering coefficient, defined 
by the average ratio of the number of links (edges) con- 
necting to its nearest neighbors of a node to the number 
of possible links between all these nearest neighbors, is 
large as similar to that in regular graphs. Large cluster- 
ing coefficient means the high frequency of "the friend 
of a friend is also his/her friend." As the SF property, 
the degree distribution follows a power-law, P(k) ~ k 1 , 
2 < 7 < 3; the fat-tail distribution consists of many 
nodes with low degrees and a few hubs with very high 
degrees. Moreover, a proposal of the universal mecha- 
nisms [4] to generate SF networks inspired to elucidate 
the topological properties. One of the advantage is that 
SF networks are optimal in minimizing both the effort for 
communication and the cost for maintaining the connec- 
tions [■>]. Intuitively, SF network is positioned between 
star or clique graph for minimizing the path length (the 
number of hops or legs) and random tree for minimizing 
the number of links within the connectivity. Another im- 
portant property is that SF networks are robust against 
random failures but vulnerable against the targeted at- 
tacks on hubs. This vulnerability called "Achilles' heel 
of the Internet" [6] frightened us. However the vulner- 
ability is a double-edged sword for information delivery 
and spreading of viruses, we expect that these properties 
will be useful for developing efficient and fault-tolerant 
networks with a defense mechanism based on the protec- 
tion of hubs. Since the SF structure is at least selected 
with self-organized manners in social and biological envi- 
ronments, the evolutional mechanisms may give insight 



to distributed network designs or social managements in 
communication or business. 

On the other hand, in contrast to abstract graphs, 
many real networks are embedded in a metric space. It 
is therefore natural to investigate the possibility of em- 
bedding SF networks in space. The related applications 
are very wide in the Internet (routers), power-grids, air- 
lines, mobile communication [7] , sensor networks [8] , and 
so on. However most of the works on SF networks were 
irrelevant to a geographical space. In this paper, focus- 
ing on the SF structure found in many real systems, we 
consider generation rules of geographical networks whose 
nodes are set on a Euclidean space and the undirected 
links between them are weighted by the Euclidean dis- 
tance. 

The organization of this paper is as follows. In section 
2, we introduce an example that shows the restriction 
of long-range links in real networks. Indeed, the decay 
of connection probability for the distance between nodes 
follows exponential or power-law. In section 3, we review 
recent studies of geographical SF network models, which 
are categorized in three classes by the generation rules. 
We refer to the analytical forms of degree distributions 
that characterize the SF structure. In section 4, we con- 
sider the relations among these models. In addition, we 
compare the properties of planarity and distance of con- 
nections. Finally, in section 5, the summary and further 
issues are briefly discussed. 



II. SPATIAL DISTRIBUTION IN REAL- WORLD 
NETWORKS 

The restriction of long-range links has been observed in 
real networks: Internet at both router and autonomous 
system (AS) levels obtained by using NETGEO tool to 
identify the geographical coordinates of 228,265 routers 
[9]. These data suggest that the distribution of link 
lengths (distance) is inversely proportional to the lengths, 
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invalidating the Waxman's exponentially decay rule [10] 
which is widely used in traffic simulations. Other ev- 
idence has been reported for the real data of Inter- 
net as AS level (7,049 nodes and 13,831 links) com- 
piled by the University of Oregon's Route Views project, 
road networks of US interstate highway (935 nodes and 
1,337 links) extracted from GIS databases, and flight- 
connections (187 nodes and 825 links) in a major airline 
[11]. It has been shown that all three networks have a 
clear bias towards shorter links to reduce the costs for 
construction and maintenance, however there exist some 
differences: the road network has only very short links 
on the order of 10km to 100km in the sharply decaying 
distribution, while the Internet and airline networks have 
much longer ones in the bimodal distribution with dis- 
tinct peaks around 2000km or less and 4000km. These 
differences may come from physical constraints in the link 
cost or the necessaries of long distant direct connections. 

As a similar example, we investigate the distributions 
of link lengths (distances of flights) in Japanese airlines 
[12]. The networks consists of 52 nodes (airports) and 
961 links (flights) in the Japan AirLines (JAL), 49 nodes 
and 909 links in the All Nippon Airlines (ANA), and 
84 nodes and 1,114 links including the international one. 
Fig. 1 shows the cumulative number of flights for the 
decreasing order of length measured by mile. We remark 
an exponential decay in the domestic flights (red and blue 
lines in the Fig. 1), while it rather follows a power-law by 
adding the international flights (green line in the Fig. 1). 
Note that the distribution of the link lengths is obtained 
by the differential of the cumulative one and that the 
decay form of exponential or power-law is invariant. 

Thus, link lengths are restricted in real systems, al- 
though the distribution may have some various forms as 
similar to the cases of degree distribution [13]. 




FIG. 1: Cumulative number of flights in Japanese airlines. 
The red, blue, and green lines correspond to domestic flights 
in JAL, ANA, and that including international flights (Inset: 
semi-log scale). The magenta and cyan lines guide the esti- 
mated functions of exponential and power-law, respectively. 
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TABLE I: Summary of geographical SF network models. The 
symbols O, A, and x denote goodness levels for each prop- 
erty. 

III. GEOGRAPHICAL SF NETWORK MODELS 

We review geographical SF network models in the 
state-of-the-art. By the generation rules of networks, 
they are categorized in three classes as shown in Table I. 
The generation rules are explained by variations in the 
balance of minimizing the number of hops between nodes 
(benefits for transfers) and the link lengths. 

In this section, we refer to the generation rules and the 
power-law behavior only in the essential forms because 
of the limited pages. The properties of planarity without 
crossing links and the link lengths will be discussed in 
the next section. 



A. SF networks with disadvantaged long-range 
links 

The modulated Barabasi-Albert (BA) model [14] and 
the geographical threshold graph [18] belong to the first 
class: SF networks with disadvantaged long-range links 
between nodes whose positions are random on a space 
[42]. They are natural extensions of the previous non- 
geographical SF network models by the competition of 
preferential linking based on the degree or weight and 
the restriction of link length (distance dependence). 

1. Modulated BA model in the Euclidean space 

Before explaining the first class, we introduce the well- 
known BA model [4] generated by the following rule: 
growth with a new node at each time and preferential 
attachment of links to nodes with large degrees (see Fig. 
2(a)). 

BA-Step 0: A network grows from an initial No nodes 
with m < No links among them. 

BA-Step 1: At every time step, a new node is intro- 
duced and is randomly connected to m previous 
nodes as follows. 

BA-Step 2: Any of these m links of the new node in- 
troduced at time t connects a previous node i with 
an attachment probability Hf A (t) which is linearly 
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proportional to the degree fej(i) of the ith node at 
time t, Uf A {t) ~ ki(t). 

The preferential attachment makes a heterogeneous net- 
work with hubs. More precisely, the degree distribution 
P(k) ~ fc~ 3 is analytically obtained by using a mean- 
field approximation [4] in the continuum approach [2], 
in which the time dependence of the degree hi of a given 
node i is calculated through the continuous real variables 
of degree and time. 

Based on a competition between the preferential at- 
tachment and the distance dependence of links, the mod- 
ulated BA model on a space with physical distance has 
been considered [14] . Note that the position of new node 
is random. The network is grown by introducing at unit 
rate randomly selected nodes on the Euclidean space (e.g. 
two-dimensional square area) , and the probability of con- 
nection is modulated according to Bij(t) ~ ki(t)l a , where 
I is the Euclidean distance between the ith (birth at time 
t) and the older zth nodes, and a is a parameter. The 
case of a = is the original BA model [4] . In the limit of 
a — > — oo, only the smallest value of I corresponding to 
the nearest node will contribute with probability 1. Sim- 
ilarly, in the limit of a — > oo, only the furthest node will 
contribute. Indeed, it has been estimated that the distri- 
bution of link lengths follows a power-law l~ s (long-range 
links are rare at S > 0), whose exponent is calculated by 
S = a + d — 1 for all values of a [14]. 

In the modulated BA model on a one-dimensional lat- 
tice (circumference), it has been numerically shown that 
[19] for —1 < a < the degree distribution is close to a 
power-law, but for a < —1 it is represented by a stretched 
exponential P(k) = aexp(— bk 1 ), where the parameters 
a, b, and 7 depend on a and m, although the SW prop- 
erty [3] is preserved at all values of a. For the transition 
from the stretched exponential to the SF behavior, the 
critical value is generalized to ct c = 1 — d in the embedded 
rf-dimensional space [14]. More systematic classification 
in a parameter space of the exponents of degree, distance, 
and fractal dimension has been also discussed [9] . 

Other related studies to the form of connection prob- 
ability ~ kf l a are the phase diagram of the cluster- 
ing properties in the a-/3 plane [15], the comparison of 
the topological properties for the special case of the con- 
nection probability proportional to the distance (a = 1, 
/3 = 0) and the inverse distance (a = — 1, j3 = 0) [16], the 
numerical investigation of the scaling for the quantities 
(degree, degree-degree correlation, clustering coefficient) 
of the network generated by the connection probability 
proportional to the degree with the exponential decay of 
the distance [17], and so on. 



2. Geographical threshold graphs 

Th geographical threshold graph [18] is a non-growing 
network model extended form the threshold SF network 
model [20] [21]. It is embedded in the d-dimensional Eu- 



clidean space with disadvantaged long-range links. We 
briefly show the analysis of degree distribution. 

Let us consider a set of nodes with the size N. We 
assume that each node i is randomly distributed with 
uniform density p in the space whose coordinates are de- 
noted by x\, X21 ■ ■ ■ , Xd, and that it is assigned with a 
weight Wi > by a density function f(w). According 
to the threshold mechanism [18], a pair of node is 
connected iff 

(tVi + Wj)h(rij) >6, (1) 

where hfcj) is a decreasing function of the distance r,j > 
between the nodes, and 9 is a constant threshold. 

If f(w) is the Dirac delta function at w* > 0, then 
the condition of connection (1) is equivalent to > 

hr x (oAr) = f f* by using the inverse function hr 1 . This 
case is the unit disk graph, as a model of mobile and 
sensor networks, in which the two nodes within the radius 
r* are connected according to the energy consumption. 
However, the degree distribution P(k) is homogeneous. 
We need more inhomogeneous weights. 

Thus, if the exponential weight distribution 

f(w) = Ae~ A ™', (2) 

and the power-law decay function h(rij) = {rij)^ 13 , (3 > 
0, are considered, then the degree is derived as a function 
of weight 

k(wi) = / f(wj)dwj / pdx\ . . . dxd ~ e Xwi , 

JO J(Mii+Mi J )/(r ij )f i >e 

(3) 

after slightly complicated calculations. The second in- 
tegral in the r.h.s of (3) is the volume of <i-dimcnsional 
hypersphere. As in Refs. [18] [21] , by using the relation of 

cumulative distributions f£ P(k)dk = f™ f(w')dw' , 
we have 

P(k) = f(w)^. (4) 

From (3) and (4), we obtain the power-law degree distri- 
bution 

P{k) ~e- 2Xw - k~ 2 . 

Note that this result is derived only if the value of (3 is 
sufficiently small, otherwise the degree distribution has a 
streched exponential decay or an exponential decay. 

On the other hand, for the power-law weight distribu- 
tion (called Parete distribution in this form) 

w* \ w J 

we similarly obtain 

k{ w )~w d/f} , p{k) - fc-a+^/rf). 
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The exponent 7 = 1 + a/3/d is a variable depends on the 
parameters a and (3. 

Furthermore, we mention a gravity model with 
h{rij) = 1/logry . In this case, the condition of con- 
nection (1) is rewritten as Wi + w j > #logrjj, and into 



WjWj 



> 



(0) 



del' 



by the variable transformations Wi = e 

\/3 dcf 



ir, 



dcf 



and (Ri 



(rij) r/9. Eq. (6) represents a physical, 



sociological, or chemical interactions with power-law dis- 
tance dependence. For the combination of (6) and the 
weight distributions f(w) in (2) and (5), we can also 
derive the more complicated forms of P(k). Thus, the 
choice of f(w) matters for the SF properties in contrast 
to an approximately constant exponent 7 w 2 in the non- 
geographical threshold graphs [21] without h{rij). 
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FIG. 2: Network generation in each model. The analytically 
obtained degree distributions for (a)-(c) follow P(k) ~ k~ 3 , 
P(k) ~ fc -7 with cutoff k c < K, and P{k) ~ k~ 1RA , -y RA ~ 3, 
respectively. 



B. SF networks embedded on lattices 

The second class is based on the SF networks embed- 
ded on regular Euclidean lattices (SFL) accounting for 
graphical properties [22] [23]. We distinguish this class 
from the first one, because the position of node is not 
randomly distributed but well-ordered on a lattice with 
a scale that gives the minimum distance. 

Let us consider a d-dimcnsional lattice of size R with 
periodic boundary conditions. The model is defined by 
the following configuration procedures (see Fig. 2(b)) on 
an assumption of power-law degree distribution. 



SFL-Step 0: To each node on the lattice, assign a ran- 
dom degree k taken from the distribution P(k) = 
Ck~ x , m < k < K, A > 2, (the normalization 
constant: C«(A- l)mr~ l for large K). 

SFL-Step 1: Select a node i at random, and connect it 
to its closest neighbors until its connectivity ki is 
realized, or until all nodes up to a distance, 



r{ki) = Ak 



i/d 



(7) 



have been explored: The connectivity quota kj of 
the target node j is already filled in saturation. 
Here A > is a constant. 

SFL-Step 2: The above process is repeated for all 
nodes. 

As in Ref. [22] , we derive the cutoff connectivity. Con- 
sider the number of links n(r) entering a node from a 
surrounding neighborhood of radius r, when the lattice 
is infinite, R — ► oo. The probability of connections be- 
tween the origin and nodes at distance r' is 



P [ k > 




P(k')dk' 



1 r' < A 

( r Y^)d(i-A) r > > A 



'(r'/A) d 

Thus, from n(r) = S d r' d - l dr' /~ /A)- P{k')dk', we ob 
tain 

d c(r/A) d 



n(r) = V d r d { - 




kP{k)dk 



P{k)dk ) 

(r/A)" J 



where Vd = Sd/d and Sd is the volume and the surface 
area of the d-dimcnsional unit sphere, respectively. The 
cutoff connectivity is then 



k c = lim n(r) = V d A d (k), 



(8) 



where (k) = J kP(k)dk denotes the average connectivity. 

If A is large enough such that k c > K, the network can 
be embedded without cutoff. Otherwise, by substituting 
(8) into (7), the cutoff connectivity k c implies a cutoff 
length 



£ = r(k c ) = (V d (k))^ d A 2 



(9) 



The embedded network displays the original (power-law) 
distribution up to length scale £ and repeats, statistically, 
at length scales larger than £. 

Whenever the lattice is finite, R < oo, the number 
of nodes is finite, N ~ R d , which imposes a maximum 
connectivity, 



(10) 



where the first approximation is obtained from 



f™P(k)dk 



_Q_ui- 

1-X K 



K 



a finite-size cutoff length is 



1/N. From (7) and (10), 
AR 1 /^" 1 ). (11) 
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These three length scales, R, £, r max , determine the 
nature of networks. If the lattice is finite, then the max- 
imum connectivity K is attained only if r max < £. Oth- 
erwise {r max > £), the cutoff k c is imposed. As long as 
mm(r ma x,C) <C R, the lattice size R imposes no serious 
restrictions. Otherwise (min(r ma2; , £) > i?), finite-size ef- 
fects bounded by i? becomes important. In this regime, 
the simulation results [22] [23] have also shown that for 
larger A the network resembles the embedding lattice be- 
cause of the rare long-range links, while the long-range 
links becomes noticeable as A decreases. 

Concurrently with the above work, Warren et al. [24] 
have proposed a similar embedding algorithm in a two- 
dimensional lattice, however the number of nodes in each 
circle is equal to the connectivity without cutoff. Thus, 
the main difference in their approaches is that a node 
can be connected to as many of its closest neighbors as 
necessary, until its target connectivity is fulfilled. 

In addition, Ref. [25] has discussed the shortest paths 
on d-dimensional lattices with the addition of an average 
of p long-range bonds (shortcuts) per site, whose length 
/ is distributed according to Pi ~ 
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is N. At time TV, a new node with degree m is added to 
the network, and if it connects to a preexisting node i, 
then the degree is updated by ki — » fcj + 1 with the equal 
probability m/N to all nodes because of the uniform ran- 
domly chosen interval. 

Thus, we have the following evolution equation, 

(in \ in 
1 - — J n(k, N) + -n(k - 1, N) + 8 k<m , 

where 8k, m is the Kronecker delta. Note that considering 
such equation for the average number of nodes with k 
links at time N is called "rate-equation approach," while 
considering the probability p(k, ti, t) that at time t a node 
i introduced at time ti has a degree k is called "master 
equation approach" [2]. 

When N is sufficient large, n(k, N) can be approxi- 
mated as NP(k). In the term of degree distribution, wc 
obtain 

-, / \ k—m 

for k 3> m {P{k) = for k < m), although it is not a 
power-law. 




FIG. 3: Growing networks with geographical attachment pref- 
erence. The blue node and links are newly added. 



C. Space-filling networks 

The third class is related to the space-filling packing in 
which a region is iteratively partitioned into subregions 
by adding new node and links. 

1. Growing small-world networks 

Let us consider the growing network with geographical 
attachment preference [2(j] as a modification of the SW 
model [3]. In this network, from an initial configuration 
with m + 1 completely connected nodes on the circum- 
ference of a circle, at each subsequent time step, a new 
node is added in an uniform-randomly chosen interval, 
and connects the new node to its m nearest neighbors 
w.r.t distance along the circumference. Fig. 3 (left) il- 
lustrates the case of m = 2. We denote n(k,N) as the 
number of nodes with degree k when the size (or time) 



(a) RAN 




(b) Deterministic AN 

FIG. 4: Apollonian Networks: (a) Random AN generated 
from an initial triangulation of square and (b) Deterministic 
AN generated from an initial triangle of back lines. The red 
and blue lines are added links at the first and second steps, 
respectively. 



2. Apollonian networks 

The growing small-world networks model [26] can be 
extended from polygonal divisions on a circle to polyhe- 
dral divisions on a sphere as shown in Fig. 3. Wc should 
remark the extended model becomes a planar graph, 
when each node on the surface is projected onto a plane 



G 



such as from a Riemannian sphere. It is nothing but a 
random Apollonian network (RAN) [27] [28], and also the 
dual version of Apollonian packing for space-filling disks 
into a sphere [29], whose hierarchical structure is related 
to the SF network formed by the minima and transi- 
tion states on the energy landscape [30]. The power-law 
degree distribution has been analytically shown in the 
RAN [27] [28]. To derive the distribution P(k), we con- 
sider the configuration procedures of RAN as follows (see 
Fig. 2(c)). 

RAN-Step 0: Set an initial triangulation with Nq 
nodes. 

RAN-Step 1: At each time step, a triangle is randomly 
chosen, and a new node is added inside the triangle. 

RAN-Step 2: The new node is connected to its three 
nodes of the chosen triangle. 

RAN-Step 3: The above processes in Steps 1 and 2 are 
repeated until reaching the required size N . 

Since the probability of connection to a node is higher 
as the number of its related triangles is larger, it is pro- 
portional to its degree as the preferential attachment. 
Thus, we have the following rate-equation 

n(k+l,N+l) = JL n (k,N)+(l - ) n(k+l,N), 

iVA V A A / 

(12) 

where the number of triangles Aa (at the grown size or 
time A) is defined as A A = 2(JV - 4) + 4 for an initial 
tetrahedron, and Aa = 2(A— 3) + l for an initial triangle, 
etc. 

In the term of P(k) n(k,N)/N, Eq. (12) can be 
rewritten as 

. „ . NkP(k) s N(k + l)P(k + l) 
(A+l P(fc+1) = ^^l+NP(k+l) ^ -p i 

iVA iVA 

By the continuous approximation, we obtain the solution 
P(k) ~ k-i RA with lRA = (N A + A)/A w 3 for large 
A. Fig. 4 (a) shows an example of RAN. 

Moreover, in the deterministic version [29] [31], ana- 
lytical forms of the power-law degree distribution P(k), 
clustering coefficient a, and the degree-degree correla- 
tion k nn {k) can be derived [29], since the calculations 
are easier in the recursive structure without randomness 
into subregions as shown in Fig 4 (b). Here, k nn (k) is de- 
fined by the the average degree of the nearest neighbors 
of nodes with degree k. It has been observed in tech- 
nological or biological networks and in social networks 
that there exists two types of correlations, namely dis- 
assortative and assortative mixings [32]. These types of 
networks tend to have connections between nodes with 
low-high degrees and with similar degrees, respectively. 
The RAN shows the disassortative mixing [29]. 

Similarly, the analytical forms in the high-dimensional 
both random [33] [34] and deterministic [35] Apollonian 



networks have been investigated by using slightly differ- 
ent techniques. They are more general space-filling mod- 
els embedded in a high-dimensional Euclidean space, al- 
though the planarity is violated. 




(a) growing SW network 




(b) growing spatial SF network 

FIG. 5: SW and SF networks generated by randomly selecting 
edges. 



3. SF networks generated by selecting edges 

Another modification of the growing SW networks [26] 
is based on random selection of edges [36] [37]. We clas- 
sify them in the relations to the partitions of interval or 
region as mentioned in 3.3.1 and a Voronoi diagram. The 
following two models give typical configurations (see Fig. 
5). 

The growing SW network generated by selecting edges 
[36] is constructed as follows. Initially, the network has 
three nodes, each with degree two. As shown in Fig. 5(a), 
at each time step, a new node is added, which is attached 
via two links to both ends of one randomly chosen link 
that has never been selected before. The growth process 
is repeated until reaching the required size N . Since all 
nodes have two candidates of links for the selection at 
each time, an exponential degree distribution has been 
analytically obtained. If the multi-selection is permitted 
for each link, it follows a power-law. The difference in 
configuration procedures for RAN is that, instead of tri- 
angulation with adding three links, two links are added 
at each step. We assume that the position of added node 
is random (but the nearest to the chosen link) on a met- 
ric space, although only the topological properties are 
discussed in the original model [36]. 

In the growing spatial SF network [37] on a two- 
dimensional space, m links whose center points are near- 
est to an added new node (as guided by the dashed lines 
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in Fig. 5 (b)) are chosen at each time step. Both end 
nodes of the nearest link(s) have an equal probability of 
connection. If a Voronoi region [38] [39] for the center 
points of links is randomly chosen for the position of new 
node in the region [43] , the selection of a link is uniformly 
random, therefore the probability of connection to each 
node is proportional to its degree. Then, we can ana- 
lyze the degree distribution. Note that any point in the 
Voronoi region is closer to the center (called generator 
point) belong in it than to any other centers. 

For the case of m = 1 as a tree, the number of node 
with degree k is evolved in the rate-equation 

(k — 1) k 
n(k,t+l) = n(k, t) + - n(k - 1, t) - —n(k, t)+S k ,i, 

(13) 

where n(k,t) denotes the number of nodes with degree 
k, and 2t is the total degree at time t. 

In the term of degree distribution P(k,t) ss n(k,t)/t 
at time t, Eq. (13) is rewritten as 

(t + l)P(k,t + l)-tP(k,t) = 

i [(k - l)P{k - 1, t) - kP(k, t)} + 6 kA . 
At the stationary value independent of time t, we have 

P(k) = \l( k - l ) p i k - 1) - kP(k)} + 6 kA . 

From the recursion formula and P(l) = 2/3, we obtain 
the solution 

P(k) = ^—^P(k - 1) = -ry, -4ti TT ~ AT 3 . 

w k + 2 K ' k(k + l)(k + 2) 

IV. RELATIONS AMONG THE MODELS 

We discuss the relations among the independently pro- 
posed models. Remember the summary of the geograph- 
ical SF network models in Table I. 

The first class is based on a combination of the pref- 
erential attachment or the threshold mechanism and the 
penalty of long-range links between nodes whose posi- 
tion is random, while the second one is on embedding 
the SF structure with a given power-law degree distri- 
bution in a lattice. Since the assigned degree to each 
node can be regarded as a fitness [2], the SFL is con- 
sidered as a special case of the fitness model [20] em- 
bedded on a lattice. In contrast, the penalty of age or 
distance dependence of each node can be regarded as a 
non-fitness in general term. If we neglect the difference 
of penalties, this explanation bridges the modulated BA 
[14] [19], SFL [22] [23], and aging models [40] with a gen- 
eralized fitness model. The crucial difference is the posi- 
tioning of nodes: one is randomly distributed on a space 
and another is well-ordered on a lattice with the mini- 
mum distance between nodes. Moreover, the weight in 



the threshold graphs [18] [21] is corresponded to a some- 
thing of fitness, however the deterministic threshold and 
the attachment mechanisms should be distinguished in 
the non-growing and growing networks. We also remark 
that, in the third class, the preferential attachment is im- 
plicitly performed, although the configuration procedures 
arc more geometric based on triangulation [27] [28] [29] or 
selecting edges [37]. In particular, the position of nodes 
in the Apollonian networks is given by the iterative sub- 
divisions (as neither random nor fixed on lattice), which 
may be related to territories for communication or supply 
management in practice. 

Next, we qualitatively compare the properties of pla- 
narity without crossing links and link lengths. We em- 
phasize that the planarity is important and natural re- 
quirement to avoid the interference of beam (or collision 
of particles) in wireless networks, airlines, layout of VLSI 
circuits, vas networks clinging to cutis, and other net- 
works on the earth's surface, etc [28]. 

In the modulated BA models and the geographical 
threshold graphs, long-range links are restricted by the 
strong constraints with decay terms, however crossing 
links may be generated. There exist longer links from 
hubs in the SFL, because such nodes have large num- 
ber of links however the positions of nodes are restricted 
on a lattice; the density of nodes is constant, therefore 
they must connect to some nodes at long distances. More 
precisely, it depends on the exponent A of power-law de- 
gree distribution as mentioned in the subsection 3.2. In 
addition, the planarity is not satisfied by the crossing be- 
tween the lattice edges and the short-cuts. On the other 
hand, RAN has both good properties of the planarity 
and averagely short links. However, in a narrow triangle 
region, long-range links are partially generated as shown 
in Fig. 4. Similarly, the SF networks generated by se- 
lecting edges may have long-range links as shown in Fig. 
5 (b): the chosen end point for connection is long away 
from the newly added node at a random position, even 
though the selected edges have the nearest centers. 



V. CONCLUSION 

In this review of geographical SF network models, we 
have categorized them in three classes by the genera- 
tion rules: disadvantaged long-range links, embedding 
on a lattice, and space-filling. We have shown that these 
models have essential mechanisms to generate power-law 
degree distributions, whose analytical forms can be de- 
rived on an assumption of the restricted link lengths as 
consistent with real data. Furthermore, the basic topo- 
logical properties of the planarity and link length have 
been discussed for each model. In particular, the geo- 
graphical threshold graphs and the RAN are attractive 
because of the tunable exponent 7 of P(k) ~ fc -7 or 
the locality related to the unit disk graphs, and the pla- 
narity of network on the heterogeneous positioning of 
nodes. However, they have drawbacks of crossing and 



8 



long-range links, respectively. To avoid long-range links, 
an improvement by the combination of RAN and Dclau- 
nay triangulation based on diagonal flipping [38] [39] is 
considering [41]. 

We have grasped several configuration procedures of 
geographical SF networks and discussed the above prop- 
erties, however these are still at the fundamental level. 
We must consider further issues, for example, 

• Quantitative investigation of the topological prop- 
erties including diameter of network, clustering co- 
efficient, degree-degree correlation, and between- 
ness centrality (related to congestion of information 
flow), etc. 



• Analysis of dynamics for the traffic and the fault- 
tolerance, especially in disaster or emergent envi- 
ronment. 



• Positioning of nodes with aggregations according 
to a population density in the evolutional and dis- 
tributed manners. 

We will progress to the next stage from the observation 
of real networks to the development of future networks. 
The distributed design and management will be usefully 
applied to many socio-technological infrastructures. 
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